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Abstract. An analysis of some properties of stellar mod- 
els as a function of metallicity Z (and helium content Y) 
is presented, with special attention to those stars with 
metallicities higher than twice or three times solar. 

It is shown that the stellar properties as a function of 
Z are mainly determined by the effects of the opacities 
at sub-solar metallicities, and by the effects of the mean 
molecular weight and stellar mass loss at higher metal- 
licities. As a result, very metal-rich stars [Z J> 0.05) ex- 
hibit properties that deviate from what is expected from 
the known characteristics at lower metallicities. In par- 
ticular, they are more luminous and hotter than those at 
Z < 0.05 due to the effect of the mean molecular weight. 
They have main sequence lifetimes much shorter (60% 
shorter at Z = 0.1 than at Z = 0.02) than those at so- 
lar metallicity due to their lower initial hydrogen content. 
Finally, the high mass loss rates at high metallicities af- 
fect significantly the population synthesis of massive stars 
in very metal-rich regions. An analysis of expected Wolf- 
Rayet and supernovae populations in such conditions is 
briefly presented. 

Key words: Stars: evolution of; Stars: general; Stars: 
Hertzsprung- Russell diagram; Galaxy (the): center of; 
Stars: Wolf-Rayet 



1. Introduction 

Observations as well as models of galactic chemical evo- 
lution suggest that the metallicity at the center of our 
galaxy, where the density of matter is high, can reach up 
to 3 to 5 times the solar metallicity (McWilliam & Rich 
1994, Simpson et al. 1995). Other objects such as elliptical 
galaxies, or maybe even quasars (Korista et al. 1996), also 
show evidences of very metal-rich stellar populations. 

The metallicity affects the evolution of stars mainly 
through its impact on the radiative opacities, the equa- 
tion of state, the nuclear reaction rates and the stellar 
mass loss rates. The impact of the metallicity on the stel- 
lar properties has been well studied for 0.001 < Z < 0.04 
on base of extensive stellar model calculations (see e.g. 



Schaller et al. 1992, hereafter Paper I, and Maeder & Conti 
1994). In contrast, they are still poorly known at metal- 
licities much higher than solar. Recent calculations have 
been performed at metallicities five times solar (Fagotto 
et al. 1994, Mowlavi et al. 1997). Mowlavi et al. (1997) 
notice that those very metal-rich models present specific 
properties which are distinct from the ones at Z < 0.04. In 
particular, they are hotter and more luminous that those 
at solar metallicity, contrary to expectations. It is the pur- 
pose of this paper to analyze in detail some of these prop- 
erties of very metal-rich stars. These could have important 
consequences on the observable properties of high metal- 
rich galaxies. 

The study is based on the extensive grids of stellar 
models computed by the Geneva group with 0.001 < Z < 
0.1. It covers the main phases of stellar evolution, exclud- 
ing the latest stages of evolution. The core helium burning 
phase of very metal-rich low-mass (M < 1.7 M Q ) stars, 
important for the study of elliptical galaxies, is not con- 
sidered in this paper as it has been the subject of existing 
papers in the literature (see, for example, Dorman, Rood 
& O'Connell, 1993, Fagotto et al. 1994). 

A general discussion on the dependence of stellar prop- 
erties on chemical composition is presented in Sect. |[ 
Section [| summarizes the model ingredients used in the 
computation of the Geneva grids. The stellar properties 
during the main sequence (MS) and the core helium burn- 
ing (CHeB) phase are analyzed in Sects. || and ||, respec- 
tively. Wolf-Rayet stars in very metal-rich environments 
are briefly discussed in Sect. [| Conclusions are drawn in 
Sect. @. 

2. Chemical composition and stellar properties 

Chemical composition is traditionally described by X, Y 
and Z, which are the mass fractions of hydrogen, helium 
and heavier elements, respectively. In the course of galactic 
evolution, the interstellar matter is progressively enriched 
in heavy elements, and Z increases through time. It is evi- 
dent that helium is also enriched, the increase of its inter- 
stellar abundance with respect to Z being conveniently de- 
scribed by a AY/AZ law such that Y = Y + AY/AZ x Z 
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(where Yq is the primordial He abundance). At low metal- 
licitics (Z <, 0.02), Y remains practically constant [an 
absolute variation 5Z of the metallicity induces an ab- 
solute variation of He of SY = AY/AZ x SZ : so that 
5Y/Y <, 1Q- 2 for Z & 10~ 3 ], and the stellar properties 
as a function of metallicity are thus expected to be de- 
termined mainly by Z. At high metallicities (Z <; 0.02), 
on the other hand, Y increases (or, alternatively, X de- 
creases) significantly with Z (Y passes from 0.30 to 0.48 
when Z increases from 0.02 to 0.1, see Table |l|). In those 
conditions, both Z and Y (and X = 1 — Y — Z) determine 
the stellar properties as a function of metallicity. 

Chemical composition influences the stellar structure 
basically through three contributions: the mean molecular 
weight /i, the radiative opacity k, and the nuclear energy 
production e nuc . The metallicity-dependent mass loss rate 
M further affects stellar evolution as a function of Z. 

2.1. The ^-effect 

The mean molecular weight acts in the hydrostatic equi- 
librium of the star. Let us imagine a thought experiment 
in which /i is arbitrarily increased over the whole star, all 
other quantities such as k or e being unaltered (though 
keeping, of course, their density and temperature depen- 
dences). The reduction in the gas pressure leads to an 
overall contraction (and heating) until hydrostatic equi- 
librium is re-established. As a result, the stellar radius de- 
creases and the surface temperature T e ff increases. The 
surface luminosity L, on the other hand, increases with \i 
since L(r) oc T 3 dT/dr [in the case of radiative transfer, 
with L(r) and T being the luminosity and temperature at 
radius r, respectively]. In summary, an increase in /i tends 
to increase both L and T e ff. 

This /^-effect is undetectable in low-metallicity stars, 
where fj, is essentially constant as a function of metallic- 
ity. At metallicities higher than about solar, however, fj, 
increases as a result of increasing Y [Z remains unim- 
portant in this respect for all realistic values of Z, since 
H ~ (2X + 0.75Y + 0.5Z)- 1 }. And indeed, it is shown in 
Sections |] and [5] that very metal-rich models are more lu- 
minous and hotter than solar-metallicity ones. It has to be 
noted that the AY/AZ law adopted in metal-rich models 
has an important influence in this respect, as it fixes Y 
and thus \x. 

2.2. The k- effect 

The radiative opacity acts in the (radiative) transport 
of energy. Let us consider again a thought experiment 
in which the opacity is arbitrarily increased throughout 
the star, all other quantities (such as [i or e nuc ) being 
unaltered. The hydrostatic structure remains unaffected 
to first approximation. In particular the radius and tem- 
perature profiles are unchanged. The energy flux, on the 
other hand, decreases with increasing k [since L(r) oc 



K ~i rp3 dT/dr f or radiative transfer], leading to a decrease 
in the surface luminosity. As the radius is unaffected to 
first order, the surface temperature decreases too. In sum- 
mary, an increase in k tends to decrease both L and T e ff. 

In general, the opacities increase with increasing 
metallicity. This is the case for the bound-free (see 
Eq. 16.107 in Cox & Giuli 1969) and free-free (see 
Eq. 16.95 in Cox & Giuli 1969) transitions. In contrast, 
electron scattering depends only on X [« e ~ 0.20(1 + X)}, 
which is about constant at Z <J 0.01 and decreases with 
Z at Z <; 0.01. As seen in Sect. these considerations ex- 
plain the stellar properties as a function of Z at sub-solar 
metallicities. 

2.3. The effect 

The nuclear energy production e nuc sustains the stellar 
luminosity. If e nuc is arbitrarily increased, the central re- 
gions of the star expand, leading to a decrease in the tem- 
peratures and densities, and to an increase in the stellar 
radius. Thus, an increase in e nuc tends to decrease both 
L and T e ff, as well as the central temperature T c and 
density p c . 

The way e nuc depends on metallicity is closely related 
to the mode of nuclear burning. When the CNO cycles are 
the main mode of H burning, e nuc increases with Z. In 
contrast, when the main mode of burning is the pp chain, 
£nuc is related to X. It is thus about independent of Z at 
Z <, 0.02, and decreases with increasing Z at Z <; 0.02. As 
shown in Sect. these properties determine the behavior 
of T c and p c in MS stars as a function of Z . 

2.4. The effect of M 

When mass loss is driven by radiation, Kudritzki et al. 
(1989) showed that M in O stars is proportional to Z 5 . 
As a result, the effects of mass loss dominate in more 
metal-rich stars. Mass loss reduces the stellar mass, mod- 
ifying thereby the post-MS stellar properties (and even 
during the MS for the most massive ones). It also turns 
the (massive) star into a Wolf-Rayet star at an earlier 
stage than it would at lower metallicity. This is analysed 
in more detail in Sect. ^. 

2.5. Summary 

At Z <, 0.02, Z is the main factor acting through the 
K-effect. The adopted AY/AZ law is thus not influential. 
At Z J> 0.02, on the other hand, Y plays a dominant role 
mainly through the /i-effect, while Z acts on the mass loss 
rate. The adopted AY/AZ law is then important. 

These combined effects of /i, k and e nuc on the stellar 
surface properties can be estimated more quantitatively 
with a semi-analytical approach using homology relations. 
This is developed in the appendix for a 3M© star. A 
more detailed analysis of stellar properties with metal- 
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Table 1. Geneva group grids of models (see text) 



Grid" 


z 


X 


Y 


M 


Interior opacities 


Low T opacities 


I 


0.001 


0.756 


0.243 


standard 


Rogers & Iglesias (1992) 


Kurucz (1991) 


III 


0.004 


0.744 


0.252 


standard 


Iglesias et al. (1992) 


Kurucz (1991) 


II 


0.008 


0.728 


0.264 


standard 


Rogers & Iglesias (1992) 


Kurucz (1991) 


I 


0.020 


0.680 


0.300 


standard 


Rogers & Iglesias (1992) 


Kurucz (1991) 


IV 


0.040 


0.620 


0.340 


standard 


Iglesias et al. (1992) 


Kurucz (1991) 


VII 


0.100 


0.420 


0.480 


standard 


Iglesias & Rogers (1996) 


Alexander & Fergusson (1994) 


V 


0.001^0.04 


M>12 


M 


2 x standard 


Iglesias et al. (1992) 


Kurucz (1991) 



w I : Schaller et al. (1992, Paper I); IV : Schaerer et al. (1993, Paper IV); 

II : Schaerer et al. (1993, Paper II); V : Meynet et al. (1994, Paper V); 

III: Charbonnel et al. (1993, Paper III); VII: Mowlavi et al. (1997, Paper VII). 



licity, however, can only be done with the help of stellar 
model calculations. Such analysis is provided in Sects. ^ 
and U using the Geneva grids of models, whose ingredients 
are first summarized in the next section. 



3. Models 

The stellar grids published by the Geneva group from 
Z = 0.001 to 0.1 (summarized in Table constitute a 
homogeneous set of models. They are computed with a 
value of a p = l/H p = 1.6 for the ratio of the mixing 
length I in convective regions to the pressure scale height 
Hp (determined in order to match the sun and the lo- 
cation of the red giant branch in the HR diagram), and 
a moderate core overshooting distance of d = 0.20 H p 
beyond the classical convective core boundary as defined 
by the Schwarzschild criterion (in order to reproduce the 
MS width in the HR diagram). We refer to Paper I (see 
Table [l] for the references of the papers) for a general de- 
scription of the physical ingredients. We mention below 
only the specific ingredients which differ from one set of 
calculations to another. 

The primordial helium abundance Yq is taken equal to 
0.24 (cf. Audouze 1987) for all models. The relative ra- 
tio of helium to metal enrichment, on the other hand, is 
AY/AZ = 3 at Z < 0.02, 2.5 at Z = 0.04, and 2.4 at 
Z = 0.1. At low metallicity, the model output is not very 
sensitive to the adopted value for AY/ AZ. At Z > 0.02, 
however, the model output is sensitive to AY / AZ (see 
Sect. |^). The values adopted at Z = 0.04 and 0.1 are con- 
sistent with the currently admitted value of ~ 2.5 (see 
Peimbert 1995 and references therein), but the uncertain- 
ties are still very high (Maeder 1998). 

The grids presented in Papers I to VII are all calcu- 
lated with the OPAL opacities. They use, however, differ- 
ent OPAL tables according to their availability at the time 
of each calculation, which are summarized in Table [j]. The 
low temperature opacities are taken from Kurucz (1991) 



in all grids, except in Paper VII which uses the Alexander 
& Ferguson (1994) tables. 

The mass loss prescription is of crucial importance 
for the evolution of the models. A general dependence of 
M oc Z 5 , as suggested by the models of Kudritzki et al. 
(1989), is adopted for all the stars except the WR stars. 
A Reimers law (1975) is adopted in all grids for low- and 
intermediate-mass stars. For massive stars, the situation 
is more delicate. In a first set of calculations (Papers I to 
IV), mass loss rates from de Jager et al. (1988, 'standard' 
mass loss rates) were used. A second set of calculations 
(Paper V) were then performed with enhanced mass loss 
rates during the pre-WR and WNL phases (by a factor 
of two over the 'standard' ones). The latter models were 
found to better reproduce many observational constraints 
concerning massive stars (Maeder & Meynet 1994), and 
we use them in this paper too. At Z — 0.1, however, the 
standard mass loss rates for massive stars are used (Pa- 
per VII). Indeed, we expect the mass loss rates driven by 
radiation pressure to be probably saturated at such high 
metallicities. The conclusions are, however, qualitatively 
not affected by this choice. 

Let us recall that at Z = 0.1, the criterion defining a 
star to be of the Wolf-Rayet type could be different from 
that at lower metallicities. However, in the absence of any 
observational counterpart of those objects at such high 
metallicities, we keep the same criterion as in Paper I, i.e. 
we consider that a star enters the WR stage when its sur- 
face hydrogen mass fraction becomes lower than 0.4 and 
its effective temperature is higher than 10000 K. These 
values are actually not crucial for the analysis of WR pop- 
ulations since the characteristic time-scale for entering the 
WR stage is very short compared to the WR lifetimes. The 
mass loss rates of these objects are then accordingly taken 
as in Paper I. 



4. The main sequence 
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Fig. 1. ZAMS locations in the HR diagram of 1, 3 and 
20 M Q star models at different metallicities as labeled in 
the figure for the 3M© star. Models at Z = 0.1 are iden- 
tified by triangles. The models are computed with OPAL 
and Alexander & Ferguson (1994) opacities 



4-1- The zero-age main sequence 

The arguments presented in Sect. enable to understand 
the location of the zero-age main sequence (ZAMS) in 
the Hertzsprung-Russell (HR) diagram as a function of 
metallicity. That location is given in Fig. [l] for stars of 1, 
3 and 20 M© at metallicities ranging from Z — 0.001 to 
0.1. 



Metallicities below ~ 0.04; At Z <, 0, .4, the surface prop- 
erties (as a function of Z) are determined by the K-effect 
(see Sect. 2.2). The surface temperature is mainly sensitive 
to the opacity in the outer layers where bound-bound and 
bound-free transitions dominate, so that all ZAMS mod- 
els get cooler with increasing metallicities. The surface 
luminosity, on the other hand, is sensitive to the interior 
opacities. In low- and intermediate-mass stars, the main 
source of opacity is provided by bound-free transitions, 
and increases with metallicity. As a result, L decreases 
with Z. In massive stars, however, electron-scattering is 
the main source of opacity. As this source is metallicity in- 
dependent, the luminosity of low-metallicity massive stars 
is insensitive to Z. 
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Fig. 2. Same as Fig. [l]but in the (\ogp c , logT c ) diagram 



Metal-rich stars: As the metallicity increases above Z ~ 
0.02, the /i-effect becomes competitive with the K-effect. In 
low- and intermediate-mass stars, the /it-effect eventually 
exceeds the K-effect above Z ~ 0.05, and a further increase 
in the metallicity results in hotter and more luminou^ 
ZAMS models (see Sect. 24). In massive stars, L already 
begins to increase as a function of Z at Z ~ 0.02 since the 
K-effect on L is negligible in those stars. 

Let us now consider the central temperature T c and 
density p c . Those quantities are, contrarily to L and 
T e ff, closely related to the nuclear energetic properties 



of the core material (see Sect. 2.3). In particular, they 
are sensitive to the mode of hydrogen burning. Usually, 
the CNO cycles are the main mode of core H burning for 
stars more massive than about 1.15 — 1.30 M© (depend- 
ing on metallicity), while the pp chain operates at lower 
masses. At Z — 0.1, however, all stars are found to burn 
their hydrogen through the CNO cycles (at least down to 
M = 0.8 M considered in our grids). 

When the CNO cycle is the main mode o f bu rning, p c 
and T c decrease with increasing Z (see Sect. 2^). This is 
clearly visible in Fig. ^ for the 3 and 20 M© stars at Z < 
0.04. At Z — 0.1, however, the higher surface luminosities 



1 These properties remained unnoticed by Fagotto et al. 
(1994), though inspection of the data available from their ta- 
bles reveal that their results indeed support our conclusions. 
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Fig. 3. Masses of convective cores 
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masses as labeled (in solar masses) next to the curves, as 
a function of metallicity 



of the models as compared to those at Z = 0.04 result in 
concomitant higher central temperatures. 

When the pp chain is the main mode of burning, on the 
other hand, the location in the (logp c , logT c ) is not very 
sensitive to Z. This is well verified for the 1 M models at 
Z<0.1 (Fig. |). 

4-2. Main sequence convective core sizes 

The mass of the convective core is mainly determined by 
the nuclear energy production, and thus by L. At Z <, 
0.04, it slightly decreases with increasing Z in low- and 
intermediate-mass stars, while it remains approximately 
constant in massive stars (see Fig. ||). At Z <; 0.04, on 
the other hand, it increases with Z at all stellar masses^]. 
We recall furthermore that at Z = 0.1 all stars with M as 
low as 0.8 M Q possess convective cores. 

4-3. Main sequence lifetimes 

The MS lifetimes tn are summarized in Fig. || for three 
different metallicities as a function of the initial stellar 
mass, and in Fig. || for several stellar masses as a func- 
tion of metallicity. They can be understood in terms of 
two factors: the initial H abundance, which determines 
the quantity of available fuel, and the luminosity of the 
star, which fixes the rate at which this fuel burns. 

2 It is interesting to note that we could have expected smaller 
convective cores in metal-rich massive stars where the main 
source of opacity is electron scattering (since this opacity is 
positively correlated with the H content). But the effect of 
higher luminosity in those stars overcomes that K-effect, and 
the core mass increases with increasing Z. 




1 10 100 



M (M a ) 

Fig. 4. Main sequence (filled circles) and core helium 
burning (open circles) lifetimes of Z=0.001 (dotted lines), 
0.02 (dashed lines) and 0.1 (solid lines) models as a func- 
tion of their initial mass 
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Fig. 5. Main sequence lifetimes of models of initial masses 
as labeled on the curves, as a function of their metallicity. 
The lifetimes are normalized for each stellar mass to its 
value at Z = 0.02 



At Z <, 0.02, tn is mainly determined by L, being 
shorter at higher luminosities. Fig. || confirms, as expected 
from L (see Fig. |l|), that tn increases with Z in low- and 
intermediate-mass stars, and is about independent of Z in 
massive stars. 

At Z > 0.02, on the other hand, the initial H abun- 
dance decreases sharply with increasing Z (the H de- 
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Fig. 6. Stellar masses at the end of the main sequence 
phase as a function of initial mass for different metallici- 
ties as labeled on the curves. Dotted lines correspond to 
models computed with the 'standard' mass loss rates for 
massive stars (Papers I to IV and VII), while thick lines 
correspond to models computed with twice that 'standard' 
mass loss rates (Paper V) 

pletion law is dictated by the adopted AY/AZ law; for 
AY/AZ = 2.4 and Y = 0.24, X drops from 0.69 at 
Z = 0.02 to 0.42 at Z = 0.1). The MS lifetimes are 
then mainly dictated by the amount of fuel available. This, 
combined with the higher luminosities at Z = 0.1, leads 
to MS lifetimes which are about 60% shorter at Z = 0.1 
than at Z — 0.02. This result is independent of the stellar 
mass for M < 40 M . Above this mass, however, the ac- 
tion of mass loss extends the MS lifetime, as can be seen 
from the 60 M© curve in Fig. |^. 

4-4- Masses at the end of the MS phase 

Mass loss determines crucially the evolution of massive 
stars, especially at high metallicities. The stellar masses 
remaining at the end of the MS phase with the 'enhanced' 
mass-loss prescription (see Sect. ^) are shown in thick lines 
in Fig. |^ (the ones with 'standard' mass-loss prescription 
is also shown in dotted lines for comparison). The most 
striking result at metallicities higher than twice solar is the 
rapid evaporation of massive stars with M > 50 M© , and 
the consequent formation of WR stars during core H burn- 
ing. At Z = 0.04, stars more massive than ^80 M© leave 
the core hydrogen burning phase with less than 5 M© . The 
same holds true at Z = 0.1 even though 'standard' mass 
loss rates are used for that metallicity. 

5. The core He-burning phase 
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Fig. 7. Tracks followed in the HR diagram during the core 
helium burning phase by 3 M models at different metal- 
licities as labeled in the figure. The points on each track 
locate the model when the core helium abundance has 
dropped to X c ( 4 He) = 0.85 

5.1. The HR diagram 

The surface properties during the CHeB phase are dic- 
tated by the same arguments developed in Sect. || and used 
to explain the MS properties (see Sect. |). At Z <, 0.04, 
both L and T e ff decrease with increasing Z, due to the 
K-effect. At Z > 0.04, on the other hand, both L and 
T e ff increase with Z, due to the /^-effect. These proper- 
ties are well verified in Figs. |7] and [s] which show the CHeB 
tracks in the HR diagram of intermediate-mass and mas- 
sive stars, respectively. 

Low- and intermediate-mass stars of low metallicity 
are known to perform a blue loop in the HR diagram and 
cross the Cepheid instability strip. The extent of that loop 
is a function of both the stellar mass (more massive stars 
having more extended loops) and the metallicity (more 
metal-rich stars having less extended blue loops). For a 
3M Q star, the blue loop is suppressed at Z <; 0.008 (see 
Fig. 0), while it is suppressed at Z > 0.02 for a 5 M© star. 
Furthermore, no blue loop is obtained for any stellar mass 
at metallicities higher than Z ~ 0.04. 

Another interesting property is the effective tempera- 
ture at which massive stars enter the CHeB phase. Low- 
metallicity massive star models ignite helium as blue su- 
pergiants, while solar metallicity models ignite helium as 
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red supergiants. This location of helium ignition as a func- 
tion of Z in the HR diagram is shown in Fig. @. As a result, 
the number ratio of blue-to-red supergiants is a decreas- 
ing function of the metallicityf]. We finally notice, again, 
the departure of Z = 0.1 models to the general rule: they 
ignite helium at a surface temperature higher than those 
of sol ar m etallicity models. This is due to the //-effect (sec 
Sect. ~ 



2J) 



5.2. Core helium burning lifetime 

The CHeB lifetime tu e is displayed in Fig. || for three dif- 
ferent metallicities as a function of the initial stellar mass, 
and in Fig. |^ as a function of metallicity for several ini- 
tial masses. It follows a pattern very similar to that of 
tn (Fig. ^|) for intermediate-mass stars, and can be un- 
derstood on grounds of the same arguments developed in 
Sect. £D| 

In massive stars, the CHcB lifetime is extended as a 
result of the lower average luminosities due to mass loss. 
As a result, the CHeB lifetime relative to its value at Z = 
0.02 increases with the initial stellar mass. 



3 Let us recall that this property is found by all existing 
stellar models, but is not supported by the observations. We 
refer to Langer & Maeder (1995) for a discussion on this issue. 
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Fig. 9. Core He burning lifetimes of models of initial 
masses as labeled on the curves, as a function of metallic- 
ity. The lifetimes are normalized for each stellar mass to 
its value at Z = 0.02 



The ratio iff/ise is displayed in Fig. 10 



Below Z ~ 

0.04, this ratio is a decreasing function of metallicity, for 
a given stellar mass, and an increasing function of stellar 
mass, for a given metallicity. Above Z ~ 0.04, this ratio 
converges to the value of tn/tHe — 5.5, irrespective of 
the stellar mass. For stars more massive than ~30M©, 
however, the high mass loss rates reduce further this ratio. 



5.3. Masses at the end of the CHeB phase 

The masses at the end of the CHeB phase of massive stars 
are displayed in Fig. |ll| for several metallicities. These are 
a decreasing function of increasing metallicity. 



6. Wolf-Rayet stars in very metal-rich 
environments 

The properties of the WR star models at metallicities 
between 0.001 and 0.04 have already been presented in 
Maeder & Meynet (1994). Here we shall focus our discus- 
sion on the expected characteristics of the WR stars in 
very metal- rich environments, i.e. at Z = 0.1. 

The minimum initial mass for a single star to become 
a WR star is a decreasing function of metallicity. It goes 
from 25 M© at Z = 0.02 to 21 M© at Z = 0.04 (see table 
1 in Maeder & Meynet 1994, 2xM case), and to approx- 
imately 17 M© at Z = 0.1. 

6.1. Mi n it ^ 60 M©: ending as helium white dwarf? 

At Z = 0.1, stars more massive than about 60 M© lit- 
erally evaporate during the main sequence phase. They 
loose so much mass that they never reach the He-burning 
phase, and will likely end their nuclear life as helium white 
dwarfs. Let us remark here that helium white dwarf stars 
are generally thought to be the end point evolution of very 
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Fig. 10. Ratio of core H to core He burning lifetimes for 
models of initial masses as labeled on the curves, as a 
function of metallicity 



low mass stars typically with initial masses between 0.08 
and 0.5 Mq. Their existence in our present universe is 
thus doubtful (unless created through binary mass trans- 
fer), since the MS lifetimes of these very low mass stars 
are greater than the present age of the universe. In con- 
trast, the above scenario (already addressed by Maeder & 
Lequeux 1982; see also Paper V), in which massive evap- 
orating stars could end their stellar life as helium white 
dwarfs, allows the formation of such objects in young star 
forming regions at very high metallicities. 

6.2. WN & WC Wolf-Ray et stars 

Stars more massive than about 50 Mq never enter the WC 
phase at Z = 0.1. Due to very high mass loss rates act- 
ing already during the MS, the H-convective core rapidly 
decreases in mass, leaving behind a great portion of the 
star with CNO-enriched material. As a consequence, the 
star goes through a very long WN phase. Moreover, these 
stars enter the core He-burning phase, if they do enter this 
phase at all, with very low masses. This implies in partic- 
ular so small mass loss rates, according to the M(M) rela- 
tion proposed by Langer (1989), that their helium burning 
core never uncovers, preventing the star to become a WC 
Wolf-Rayet star. 

6.3. Wolf-Rayet versus O type stars 

The number ratio of WR to O type stars expected in a 
constant star formation rate region at Z = 0.1 is esti- 
mated from these models to be around 0.45. To compute 
this value, only the single star evolutionary channel was 
considered, we chose a Salpeter initial mass function (i.e. 
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Fig. 11. Stellar masses at the end of the core helium burn- 
ing phase as a function of initial mass for different metal- 
licities as labeled on the curves. Only models computed 
with twice the 'standard' mass loss rates for massive stars 
are shown, except at Z = 0.1 for which the mass loss rates 
are 'standard' 

with a slope equal to 1.35), and the upper mass limit was 
taken equal to 60 Mq. 

From the above result, one expects to find at Z = 0.1 
about 1 WR star per each 2 O-type stars detected! This is 
a very large proportion compared to regions with solar and 
sub-solar metallicity (in the solar neighborhood, the ob- 
served number ratio WR/O is 0.1). This high proportion 
results mainly from the very high mass losses experienced 
by these metal-rich stars, which, as seen above, decrease 
the minimum initial mass for WR star formation, make 
the star enter the WR phase at an earlier stage of its 
evolution, and increase the core helium burning lifetime. 
Moreover, the duration of the O-type phase of metal- rich 
stars is considerably reduced as compared to that at lower 
metallicities. For instance when Z increases from 0.04 to 
0.1, the lifetime of the O-type phase decreases by about a 
factor 2. This is essentially a consequence of the shorter 
duration of the H-burning phase at Z = 0.1 (see Fig. ||). 

6.4- Supernovae 

Let us now consider the final evolution of massive stars 
and estimate the relative fraction of supernovae with WR 
progenitors using some simplifying assumptions (see also 
Meynet et al., 1998). There is at present some uncer- 
tainties concerning the fate of WR stars, and to which 
type of supernovae they can give birth. As a first ap- 
proximation, we make the hypothesis that all WR stars 
do lead to supernovae events (see a discussion in Maeder 
& Lequeux 1982). Moreover, since the H-rich envelope is 
ejected through the wind during the WNL phase, we fur- 
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ther assume that WNE and WC stars result in super- 
novae of type I (whether lb, Ic, or another yet unknown 
subclassification). With these assumptions, and not con- 
sidering type la supernovae which have low-mass progeni- 
tors, one expects that metal-rich regions be characterized 
by a higher fraction of type I(b/c) supernovae relative to 
the total number of supernovae (Maeder 1992). As a nu- 
merical example, using the same IMF as in Sect. 6.3 and 



making the hypothesis of constant star formation rate, 
one obtains that, at Z = 0.1, the fraction of supernovae 
with WNE-WC progenitors amounts to 32% of the total 
number of supernovae having massive star progenitors. In 
other words, 1 supernova event out of 3 should be of type 
I(b/c) in very metal-rich regions. For comparison, this pro- 
portion decreases to 1 WR explosion per 5 — 7 supernovae 
at solar metallicity. 

7. Conclusions 

This paper analyzes some properties of stellar models as 
a function of the metallicity, with particular attention to 
those at very high metallicities (above twice or three times 
solar). The stellar properties as a function of Z are shown 
to result from the action of four ingredients^: the mean 
molecular weight, an increase of which leads to more lu- 
minous and hotter stars (the /x-effect); the opacity, an in- 
crease of which leads to less luminous and cooler stars 
(the K-effect); the nuclear energy production (but to a 
lesser extent), an increase of which leads to less luminous 
and cooler stars with lower T c and p c (the £ nuc -effect; and 
mass loss, which acts indirectly on the stellar structure by 
reducing the stellar mass (dramatically in the most mas- 
sive metal-rich stars). The /{-effect plays the main role 
at Z <J 0.04, while the /x-effect becomes preponderant at 
Z > 0.02. In the later case, the AY/ AZ law is important. 

Based on these considerations and on the Geneva's 
grids of stellar models, it is shown in particular that very 
metal-rich stars exhibit properties which cannot be ex- 
trapolated from model characteristics at lower metallici- 
ties. In particular, 

- very metal-rich stars (such as at Z = 0.1) are more lu- 
minous and hotter than those at Z < 0.05. This property 
contrasts with the usually known trend at lower metallic- 
ities (whereby L and T e f / decrease with increasing Z due 
to the K-effect); 

- very metal-rich stars have much lower MS lifetimes than 
those at Z < 0.05. For example the MS lifetime at Z = 0.1 
is 2.5 times shorter than that at solar metallicity, indepen- 
dent of their initial mass (at least up to 40 M Q ); 

- stellar population synthesis in metal-rich environments 
is significantly affected by the increased mass loss rates at 
high Z. In particular, a) massive stars with Z = 0.1 enter 

4 We neglect in this paper the effects of mild mixing processes 
due for example to rotation, diffusion, gravitational settling, 
etc... which may modify the internal composition profiles and 
some of the values given here 



the WNL Wolf-Rayet phase very early during their evo- 
lution, and those more massive than about 50 M Q prob- 
ably never enter the WC phase; b) the most massive 
(M > 60 Mq at Z = 0.1) high-mctallicity stars loose the 
majority of their mass during their MS phase, and prob- 
ably end their life as He white dwarfs; and c) the number 
of WR stars over O-type stars increases with metallicity, 
and metal-rich regions are expected to be characterized 
by a high fraction of type I (b/c/other?) supernovae rel- 
ative to the total number of supernovae having massive 
star progenitors (up to 1 out of 3 core collapse supernovae 
events at Z = 0.1). 
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Appendix: The ZAMS location in the HR diagram 
as a function of metallicity: a semi-analytical ap- 
proach 

The luminosity and effective temperature on the zero age 
main sequence (ZAMS) are given by the homology rela- 
tions (see, for example, Cox and Giuli 1969, p 696) 



L oc 



and 



l eff 



-1.02 
-1.08 



-0.03 
-0.10 



(•1) 



(•2) 



where p, and M are the mean molecular weight and the 
stellar mass, respectively, eo the temperature and density 
independent coefficient in the relation for the nuclear en- 
ergy production e nuc — e p x T u (see Cox and Giuli 1969, 
p 692), and kq the opacity coefficient in the Kramers law 
k = k,q p T~ 3,5 . The numbers on the first line of the su- 
perscripts in Eqs. [l] and[^ apply when the nuclear energy 
production results from the pp chain, while the second line 
applies when the CNO cycles provide the main source of 
nuclear energy. 

A variation AZ in the metallicity of a star of mass M 
affects its position in the HR diagram by 



A log L = 



7.3 
7.8 



A log// 



1.02 
1.08 



A log k 



0.02 
0.08 



Alogeo (.3) 



AlogT e// = {^}Alo g/ ,-{^}Al og 



K 



0.03 
0.10 



Alogeo (-4) 
The variation of \i is easily related to that of Z: 



Alog M =i| M A^ 

where we have used the relations 

Y = 0.24 + 2.4Z 

and 

X + Y + Z = 1. 
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Fig. .1. (a): Density profile as a function of temperature 
in a ZAMS model of a 3M Q stars at Z = 0.1; (b): Opacity 
profiles, as a function of temperature, at three metallicities 
as labeled on the curves, and calculated from the tables 
of Iglesias & Rogers (1996) at the densities displayed in 
(a); (c): Same as (b), but for the derivative dlogre/dZ 
at the labeled metallicities. The derivative is calculated 
numerically from (logK Z+0 01 - log k z ~ 001 )/0.02 



The variation of A log kq as a factor of A Z is more diffi- 
cult to evaluate. We can get an estimation of A log kq/AZ 
by comparing the opacities given by the Rogers and Igle- 
sias (1992) tables at different metallicities. The opacities 
at Z = 0.02, 0.06 and 0.1 (and at X,Y given by Eqs. [| 
and ^) are shown in Fig. y](b) for temperature and den- 
sity conditions representative of the interior of a 3M 0) 
Z = 0.1 ZAMS model [Fig. [|fa); we assume, as a first ap- 
proximation, that the same p — T profile applies at those 
three metallicities]. If the opacity profile were following a 
Kramers law, then the derivative d log n/dZ would be con- 
stant and equal to that of the Kramers opacity coefficient 
dlogKo/dZ (p and T being held constant). Inspection of 
Fig. [l](c) shows that this is obviously not the case. The 
derivative is seen to vary from 4 to 15 at Z — 0.02 and 
from 1.5 to 3.3 at Z = 0.09. However in this order of 
magnitude approach, we consider mean values equal 



A log kq ~ 10AZ 
AlogKo ^ 2.5 AZ 



Z = 0.02 
Z = 0.1 



(.8) 
(.9) 



The variation of A log eo as a function of heavy element 
abundances AZ can be estimated knowing that the CN 
cycle is the main energy producer in the 3 M Q star (see e.g. 
Mowlavi 1995, Fig. A.3). The two reactions 12 C (p , 7) 13 N 
and 13 C (p , 7) 14 N contribute for more than 70% of the 
total nuclear energy at the center of the star. We can thus 
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write that eq <x Z x X, which leads to, using Eqs. [6| and 



A,„ g<0 =(l-|f)J_ AZ ( . 10) 

Equations [3] and [4| become with Eqs. [^, [8|, and 
If], and taking X = 0.70, /i = 0.617 for Z = 0.02 and 

(.11) 



X = 0.42, n = 0.80 for Z = 0.1, 

Alogi f 0.41-0.50-0.07 = -0.16 Z = 0.02 



AlogZ \ 2.81-0.62-0.02 = 2.17 Z = 0.10 

AlogT e// _ f 0.09-0.16-0.09 = -0.16 Z = 0.02 . 
AlogZ ~~ \ 0.79-0.20-0.02 = 0.57 Z = 0.10 ^ 



We note from Eq. .11 the small effect of the nuclear pro- 
duction rate on the surface luminosity. 

These relations should not be considered very accurate 
given the simplifications made in their derivation (mainly 
because of the departure of the opacity profile from the 
Kramers law and the resulting uncertainty in Eqs. ^ and 
[9|). They however enable to understand qualitatively the 
ZAMS location of a 3M Q star in the HR diagram as a 
function of the metallicity, and support the following con- 
clusions found in Fig. |l|: 

- The surface luminosity at ZAMS is a decreasing function 
of Z at Z — 0.02 (and lower metallicities) and an increas- 
ing function of Z at Z = 0.1. The metallicity at which the 
inversion operates can only be determined from complete 
stellar model calculations and is found to be Z ~ 0.05. 

- The absolute shift |AlogL| resulting from an increase 
AlogZ is much greater at Z = 0.1 than at 0.02 (by a 



factor of about 14 according to Eq. .11). This result is 
qualitatively well verified by the ZAMS model calculations 
displayed in Fig. [j} 

- The effective temperature is a decreasing function of Z 
for Z <j 0.05, but it increases with metallicity for Z ;> 0.05. 



